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Preface : 

This  report  develops  an  analysis  of  the  transient  response 
of  a  gun  barrel  subject  to  repeated  firing.   A  finite-element 
approach  is  used.   A  piecewise  cubic  approximation  of  the  de- 
flection ensures  continuity  of  deflections  and  slopes  at  the 
net  points.   Initially,  linear  and  angular  displacements  and 
velocities  are  prescribed.   This  arbitrariness  in  the  initial 
conditions  allows  for  residual  motions  from  the  preceding  shot. 
A  digital  computer  program  based  upon  the  analysis  has  been  de- 
veloped.  The  program  may  be  used  for  the  period  during  which 
the  projectile  leaves  the  muzzle  by  setting  terms  pertaining  to 
the  projectile  and  the  gas  pressure  equal  to  zero.   If  desired, 
the  program  may  be  used  to  determine  the  static  deflection  of 
the  gun  barrel.   Horizontal  transverse  vibrations  of  the  barrel 
may  be  treated  by  discarding  the  acceleration  due  to  gravity. 
The  projectile  is  treated  as  a  point  mass.   Consequently,  pro- 
jectile balloting  is  not  included  in  the  theory.   The  effects 
of  a  gun-barrel  tuning  mass,  eccentricity  of  the  breech, 
eccentricity  of  the  recoil  mechanism,  stiffness  and  damping  of 
the  supports,  and  prescribed  motion  of  the  gun  foundation,  as 
well  as  Bourdon  pressure  and  projectile-barrel  friction  are 
included  in  the  model.   In  the  computer  program,  the  time  in- 
tegration routine  employs  Newmark's  beta  method.   The  computer 
program  listing  is  not  included  here.   However,  information  on 
the  computer  program  is  available  by  writing  to 


Professor  A.  P.  Boresi 

Department  of  Theoretical  and  Applied  Mechanics 

University  of  Illinois 

Urbana,  Illinois   61801 


NOTATIONS 

t    =   time.   Dot  denotes  time  derivative 

x   =   axial  distance  along  the  barrel  from  the  breech 
(Fig.  1) 

y   =   deflection  of  the  barrel;  y  =  3y/9t 

8    =   notation  of  cross  section  of  barrel.   If  shear 
deformation  is  discarded,  9  =  9y/3x  =  y' 

L  =  length  of  barrel 

a  =  radius  of  the  bore 

n  =  coordinate  of  the  tuning  mass 

mT  =  mass  of  the  tuning  mass  (Fig.  2) 

I_,   =   moment  of  inertia  of  the  tuning  mass  about  a  diametral 
axis  through  its  center  of  mass  (Fig.  2) 

m   =   mass  of  the  projectile 

M   =   mass  of  the  breech  block 

M   =   mass  of  the  barrel 

m        =  mass  of  the  part  of  the  barrel  in  the  range  greater 
than  x  (Eq.  26) 

T   =   moment  of  inertia  of  the  breech  block  about  a  trans- 
verse axis  through  its  center  mass 

g   =   acceleration  of  gravity 

E    =   Young's  modulus 

e,d  =   lengths  defining  the  eccentricity  of  the  breech 
(Fig.  1} 

a    =   angle  that  the  barrel  forms  with  the  horizontal 
CFig.  1) 

a„   =   spring  constant  for  the  recoil  mechanism  (assumed 
constant 

u„   =   viscous  damping  coefficient  for  the  recoil  mechanism 

v    =   v(t)  =  speed  of  the  projectile  relative  to  the  breech 

£    =   5(t)  =  displacement  of  the  projective  relative  to  the 
breech.   v  =  t, 


P   =   P(t)  =  gas  pressure  behind  the  projectile 

F    =   F(t)  =  axial  frictional  force  of  the  projectile  on 
the  barrel 

u   =   u(t)  =  recoil  displacement  of  the  barrel  (positive 
backwards,  Fig.  1) 

P   =  P(t)  =  tension  in  the  barrel  at  section  x  due  to 
inertia  (Eq.  26) 


0 


P 


=   mass  density  of  barrel 


x,  ,  x„ ,  ....  x  ,,  =  net  points.   x.,  =0,  x  ,,  =  L 
1'   2        n+1       "         1     '   n+1 

A.  =   x  .  , ,  -  x  .  .   This  is  the  length  of  the  j-th  interval 
J       J+l     J 

yj  =  yCxj'  tJ 

a-  =   constant  for  translat ional  spring  support  at  point  x. 

a'.  =   constant  for  rotational  spring  support  at  section  x. 

u.  =   constant  for  translation  damper  at  point  x. 

u  I  =   constant  for  rotational  damper  at  point  x. 


-      value  of  P  at  x  =  x-j  .   Note  that  P  .  depends  on  the 
location  %    of  the  projectile      ^ 


9.   =  8.(t)  =  value  of  6  at  x  =  x  . 

J  J  J 

r    =  shear  coefficient  (Eq.  8) 

B Jn  =  coefficients  for  the  Bourdon  effect  (Eqs.  24  and  25) 

AJ .  =  coefficients  for  the  axial  inertia  effect  (Eqs.  28 

a*  and  31) 

x_,   =  net  point  ahead  of  the  tuning  mass.   x_,  <  n  <  x_ 

1+1  (Fig.  2)  l  i+i 

e',e"  =  eccentricities  of  breech  spring  and  damper  (Fig.  1) 

p    =  p(x)  =  mass  of  barrel  per  unit  length 

n    =  number  of  segments;  n+1  =  number  of  net  points 


1.   INTRODUCTION 

A  finite-element  method  is  used  to  analyze  the  transient 

response  of  a  gun  barrel  due  to  repeated  firing.   The  dynamic 

behavior  of  the  barrel  may  be  computed  for  the  period  in  which 

the  projectile  is  in  the  barrel;  of  particular  importance  are 

the  slope  and  the  lateral  velocity  at  the  muzzle.   By  setting 

to  zero  terms  pertaining  to  the  projectile  and  the  gas 

pressure,  one  may  follow  the  gun  response  also  for  the  period 

after  the  projectile  leaves  the  muzzle.   The  barrel  is  modeled 

as  a  tapered  elastic  beam  with  finite  degrees  of  freedom  by 

subdividing  it  into  n  intervals  (finite  elements)  with  n+1 

net  points  (nodes).   The  generalized  coordinates  are  the 

deflections  y.  and  the  rotation  9.  of  the  barrel  cross  section 
J  J 

at  the  nodes.   Hence,  there  are  2(n+l)  generalized  coordinates. 
Two  additional  coordinates  that  enter  the  equations  are  the 
axial  recoil  displacement  u  of  the  barrel  and  the  axial  dis- 
placement Z,   of  the  projectile  relative  to  the  breech.   However, 
these  coordinates  do  not  enlarge  the  system  of  differential 
equations  if  the  axial  movements  of  the  barrel  and  the  projec- 
tile are  regarded  as  known  functions  of  time. 

The  breech  is  treated  as  a  rigid  mass  with  its  center 
offset  from  the  axis  of  the  barrel.   Viscous  dampers  (dashpots) 
may  be  located  at  all  nodes  or  they  may  be  provided  at  only  a 
few  points,  since  in  the  computer  program,  a  damper  can  be 
eliminated  by  setting  its  coefficient  equal  to  zero.   Also  the 
offset  of  the  breech  may  be  set  equal  to  zero  if  the  center 
of  mass  of  the  breech  lies  on  the  axis  of  the  barrel.   The 


structure  supporting  the  barrel  and  the  breech  is  modeled  as 
a  set  of  independent  translational  and  rotational  springs  and 
dashpots  at  the  net  points.   At  nodes  where  there  is  no  such 
support,  the  spring  constants  and/or  dashpot  constants  are 
set  equal  to  zero.   A  tuning  mass  may  be  located  at  an  arbitrary 
point  along  the  barrel. 

The  foundation  upon  which  the  gun  barrel  and  its  supports 
rest  is  regarded  as  a  rigid  body  that  moves  in  a  prescribed 
way.   In  general,  it  has  six  degrees  of  freedom.   However,  in 
a  linear  theory,  oscillations  of  the  barrel  in  a  vertical  plane 
are  uncoupled  from  transverse  horizontal  displacements.   Also, 
the  rotation  of  the  barrel  about  its  axis  is  uncoupled.   Con- 
sequently, attention  may  be  restricted  to  a  plane  motion  of  the 
barrel.   Oscillations  of  the  barrel  in  a  fixed  vertical  plane 
are  considered.   Horizontal  lateral  oscillations  may  be  super- 
imposed.  The  program  for  horizontal  lateral  oscillations  is 
the  same  as  for  oscillations  in  a  vertical  plane,  except  that 
effects  of  gravity  are  eliminated. 

The  differential  equations  of  motion  may  be  obtained  from 
Hamilton's  principle,  as  the  Lagrangian  equations  for  a  noncon- 
servative  system.   In  terms  of  the  virtual  work  5W  and  the 
kinetic  energy  T,  Hamilton's  principle  states  that 

h 

J  (6T  +  6W)dt  =  0  (1) 

4 

where  t<  denotes  time.   However,  because  certain  contributions 
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to  the  differential  equations  of  motion  can  be  written  down 
directly  in  terms  of  well-known  finite  element  procedures,  it 
is  expedient  to  do  so.   Hence,  in  this  report  Hamilton's 
principle  is  used  mainly  as  a  guide  to  verify  the  form  of 
certain  unusual  effects  such  as  the  Bourdon  effect ,  the  axial 
inertia  effect  and  other  pressure  effects  that  are  inherent  in 
gun  systems. 

Contributions  to  6W  come  from  several  sources.   One  source 
is  the  action  of  gravity  on  the  breech,  the  barrel,  the 
projectile,  and  the  tuning  mass.   Another  source  is  the  axial 
frictional  force  of  the  projectile  on  the  barrel.   It  introduces 
a  nonconservat ive  effect.   Another  contribution  comes  from  the 
so-called  "Bourdon  effect",  which  arises  because,  in  the  bent 
barrel,  the  area  of  the  bore  above  the  neutral  plane,  on  which 
gas  pressure  acts,  is  slightly  greater  than  the  area  below  the 
neutral  plane.   Still  another  effect  contributing  to  6W  is  the 
axial  inertia  of  the  recoiling  barrel,  which  initially  exerts 
a  straightening  influence.   The  contribution  of  the  strain 
energy  of  the  barrel  due  to  flexural  stresses,  etc.  may  be 
introduced  as  a  negative  contribution  to  5W.   This  latter  effect 
essentially  results  in  a  finite  element  model  of  the  gun  barrel 
in  the  final  differential  equation  form  for  the  gun  system. 
Accordingly,  a  finite  element  beam  model  is  used  to  incorporate 
this  effect  in  the  differential  equations,  rather  than  derive 
this  effect  from  Eq.  (1)  directly.   Finally,  to  account  for  the 
effects  of  the  supporting  structure,  the  support  structure  is 
regarded  as  sets  of  linear  rotational  and  translat ional  decoupled 
springs  and  dashpots  located  at  the  nodes  of  the  finite  element 
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model  of  the  barrel.   These  elements  are  assumed  to  be  attached 

to  the  gun  foundation  which  may  undergo  prescribed  motions. 

The  spring  constants  and  the  damping  constants  for  the  trans- 

lational  and  rotational  springs  and  dashpots  at  node  j  are 

respectively  (a.,  a'.)  and  (u  . ,  u').   At  a  node  where  there  is 
J    J         3         3 

actually  no  supporting  spring  or  dashpot ,  the  spring  constants 
and/or  dashpot  constants  are  set  equal  to  zero  in  the  computer 
program. 

In  general,  for  a  finite  element  model  of  the  barrel,  the 
net  virtual  work  may  be  expressed  as 

n+1  n+1 

6W  =  ^    Q.5y.  +  V^   Q'.59.  +  R6u  +  S5£       (2) 
L—d  3       J    L-u         3       3 

j  =  l  j  =  l 

where  n  =  the  number  of  elements  (n+1  =  number  of  nodes)  and 
where  Q.,  Q '.  ,  R",  and  S  are  the  components  of  generalized  force; 
i.e.,  the  coefficients  of  the  generalized  coordinates  u,  E,  ,    y., 
6.  in  the  expression  for  5W. 

The  kinetic  energy  T  is  the  sum  of  the  kinetic  energies  of 
the  breech,  the  projectile,  the  barrel,  and  the  tuning  mass. 
It  is  a  function  of  the  generalized  coordinates  u,  £,  y.,  8., 
where  dots  denote  derivatives  with  respect  to  time  t.   Conse- 
quently, the  differential  equations  (the  Lagrangian  equations) 
of  motion  are,  by  Hamilton's  principle, 


_d_  3T  3T  QJ_1T     JT 

dt   •  "  3y  "  V  dt  39   "  39   "  ^j  • 

3  3 

3    =  1,  n+1.  (3) 
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d     3T         3T   _    D  d     3T         3T  _  ,  „  ,         .  ,  , 

-t—  —   -  r~  =   R,  -rr  — r   -   t?  =  s  .  (3)cont'd 

dt    ~  •  3u  dt    r,  i         9£  v 

dU  d  c, 


where  Q.,  Q'.  ,  R,  S  are  the  components  of  generalized  force  (Eq.  2) 
The  Lagrange  equations  for  u  and  E,    are  irrelevant,  if  we  know 
u(t)  and  £(t)  as  functions  of  time  t.   It  is  to  be  noted  that 
3T/3y   =  3T/36   =  0,  since  T  =  T(u,  E,  ,    y    8  ). 

By  the  usual  procedures  of  finite  element  methods,  Eqs .  (3) 
may  be  written  in  the  form 

Mx  +  Cx  +  Kx  =  f(t)  (4) 

where  M  =  M(t),  C  =  C(t)  and  K  =  K(t),  are  the  mass,  damping 
and  stiffness  matrices  for  the  gun  system,  f(t)  is  the  driving 
force  vector  and  x  is  the  vector  of  generalized  coordinates, 

x  =  (y^  81,  y2,  62,  ...,  yn+1,  9n+1}T  .        (5) 

The  Lagrange  equations  (Eqs.  3  or  4)  are  linear  second  order 
differential  equations  for  the  unknown  functions  y.(t)  and  9.(t). 
They  are  adaptable  to  various  integration  methods,  such  as  the 
Runge-Kutta  method,  Newmark's  beta  method,  Houbolt's  method, 
etc.   Discontinuities  occur  at  the  instantaneous  location  of 
the  projectile,  since,  for  example,  the  gas  pressure  behind  the 
projectile,  causing  the  Bourdon  effect,  terminates  there. 
However,  discontinuities  do  not  obstruct  the  finite  element 
method.   This  is  an  advantage  over  the  direct  approach  via  the 
differential  equations  of  beams,  since  with  the  latter  formula- 
tion, the  action  of  the  projectile  is  described  with  a  Dirac 
delta  function  or  a  Heaviside  step  function.   These  difficulties 
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are  avoided  by  the  finite-element  treatment.   A  computer  program 
based  on  Newmark's  beta  method  has  been  devised  for  solving  the 
differential  equations  for  arbitrary  initial  values  y.(0),  8.(0) 
and  y.(0),  9.(0).   Vertical  bending  and  horizontal  bending  of 
the  gun  barrel  are  decoupled  in  linear  theory.   Both  of  these 
types  of  deformation  are  covered  by  the  program.   For  horizontal 
motion,  there  is  no  gravitational  effect. 

The  program  provides  the  history  of  the  motion  of  the  barrel 
from  the  initial  instant  at  which  the  projectile  is  fired,  until 
the  projectile  leaves  the  muzzle.   In  addition,  it  may  be  used 
for  the  period  after  the  projectile  has  left  the  muzzle.   To 
account  for  this  latter  period,  all  terms  pertaining  to  the 
projectile  (e.g.,  gas  pressure,  projectile  mass,  friction)  are 
dropped . 

To  implement  the  program,  we  must  know  the  recoil  displace- 
ment u(t),  the  projectile  displacement  |(t)  along  the  barrel, 
the  base  pressure  p(t)  on  the  projectile,  and  the  frictional 
force  F(t),  between  the  barrel  and  the  projectile.   If  5 ( t )  and 
p(t)  are  known,  F(t)  is  determined  by  Newton's  law  applied  to 
the  axial  motion  of  the  projectile.   In  applying  the  program 

to  an  actual  gun  system,  the  spring  constants  (a.,  a'.)  and  the 
°  *    6  3         3 

damping  constants  (u.,  u '.  )  of  the  supporting  structure  must  be 

estimated.   Because  of  the  complicated  nature  of  the  structure 

and  the  peculiar  damping  devices  in  use,  this  may  be  difficult. 

Of  particular  significance  for  the  accuracy  of  shooting  are 

the  muzzle  displacement,  y  , , .  muzzle  rotation  9  L,  and  the 

'  Jn+1  n+1 

muzzle  velocity  y  ,,  at  the  instant  that  the  projectile  leaves 

n+1  ° 
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the  barrel.   The  program  permits  a  study  of  these  quantities 
corresponding  to  various  initial  motions  and  deflections  of 
the  barrel  at  the  instant  of  firing. 

2.   DESCRIPTION  OF  THE  MATHEMATICAL  MODEL 

Figure  1  represents  a  gun  barrel  with  mass  M  and  elevation 
angle  a.   Axes  (x,  y)  move  axially  with  the  barrel  as  it  recoils, 
but  they  undergo  no  lateral  displacement.   The  projectile  is 
regarded  as  a  point  mass  m.   Therefore,  balloting  of  the  projec- 
tile is  not  considered.   A  short  time  t  after  the  projectile 
is  fired,  it  lies  at  the  barrel  point  x  =  ^(t),  as  shown  in 
Fig.  1.   The  axial  velocity  of  the  projectile  relative  to  the 
barrel  is  v(t)  =  E, .      The  axial  frictional  force  of  the  projectile 
on  the  barrel  is  FCt).   The  gas  pressure  driving  the  barrel  is 
pCt).   There  is  a  tuning  mass  m_  with  constant  coordinate  n. 
Its  center  of  mass  is  considered  to  lie  on  the  axis  of  the 
barrel.   Its  moment  of  inertia  about  a  transverse  axis  through 
its  center  of  mass  is  I    The  breech  is  regarded  as  a  rigid 
block  with  mass  M  and  moment  of  inertia  I  about  a  transverse 
axis  through  its  center  of  mass.   The  location  of  the  center  of 
mass  is  defined  by  two  lengths,  e  and  d  shown  in  Fig.  1.   The 
axial  recoil  displacement  is  u(t).   The  breech  recoil  spring 
and  recoil  dashpot  are  assumed  to  be  linear.   Their  constants 
are  0*r»  >  yR)  and  their  eccentricities  are  (e',  eM),  as  shown 
in  Fig.  1. 

If  there  is  initial  bending  of  the  barrel  due  to  weight, 
unsymmetrical  thermal  gradients,  or  manufacturing  imperfections, 
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the  bending  is  aggravated  by  inertial  reaction  between  the  barrel 
and  the  projectile.   Axial  friction  F(t)  also  excites  motion  of 
the  barrel.   Another  disturbance  comes  from  the  so-called 
"Bourdon"  effect,  which  arises  because,  in  the  bent  barrel, 
the  area  of  the  bore  above  the  neutral  plane,  on  which  the  gas 
pressure  acts,  is  slightly  greater  than  the  area  below  the 
neutral  plane.   Still  another  effect  contributing  to  the  bending 
of  the  barrel  is  axial  inertia  of  the  recoiling  barrel  and  breech 

3.   CONTRIBUTIONS  OF  VARIOUS  PARTS  OF  THE  SYSTEM  AND  OF  VARIOUS 
EFFECTS  TO  THE  EQUATIONS  OF  MOTION 

3.1.   Strain  Energy  of  the  Barrel.   The  barrel  may  be 

regarded  as  a  tapered  elastic  beam  of  length  L.   If  bending 

strain  energy  only  is  included,  the  strain  energy  per  unit  length 

is  proportional  to  the  square  of  the  curvature.   If  the  effect 

of  transverse  shear  deformation  is  included,  an  appropriate 

term  for  the  shear  energy  must  be  added.   In  this  report  these 

energies  are  taken,  respectively,  as  (1) 


UD   ,.    =  i  /    EI(y")2dx  (6) 

Bending   2  /      w  J  ' 

o 


L 
Shear     2 


Uc,a_    =  \    /   GA(y'-9)2dx  (7) 


o 

where  primes  denote  derivatives  with  respect  to  x,  y  =  the 
lateral  displacement  of  the  barrel,  0  =  the  rotation  of  the 
barrel  section  (Note  that  if  shear  deformation  is  not  included 
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the  rotation  9  =  y',  the  slope  of  the  axis  of  the  barrel), 

EI  =  the  bending  rigidity  and  GA  =  the  shear  rigidity  of  the 

beam.   By  Hamilton's  principle,  these  energies  contribute 

terms  to  the  components  of  generalized  forces  (Eq.  3)  or 

equivalently  to  the  stiffness  matrix  K  (Eq.  4). 

In  the  usual  method  of  finite  elements  (2),  the  barrel  is 

divided  into  n  intervals  (elements)  of  length  A,,  A0,  ...,  A., 

1    «         3 

.  .  .  ,  A  ,  by  points  x,  ,  x0 ,  . . .  x   , ,  where  x   ,  is  the 
n  ±         _        n  +  i.  n"f".L 

coordinate  of  the  muzzle  (Fig.  1).   These  points  need  not  be 
equally  spaced.   Then  with  y  taken  as  a  cubic  polynomial  in  x 
in  each  interval  of  length  A . ,  and  9  taken  as  a  corresponding 
quadratic  polynomial  in  x,  the  element  stiffness  matrix  for  the 
jth  element  is  (see  Appendix  I),  for  j  =  1,  2,  ...,  n, 


k. 
-3 


E.I. 
J  3 


A  .  (1+r) . 

3  3 


12 


6A 


-12 


-12 


-6A 


12 


6A 


6A   A  (4+r)   -6A    A  (2-r) 


-6A 


6A   A  (2-r)   -6A    A  (4+r) 


(8) 


where  r  =  12EI/(AGL  )  and  A.  =  x  .  , ,  -  x..   Appropriate  additions 

3  J  +  l     J 

of  the  matrices  k .  are  required  to  incorporate  the  element 
stiffnesses  into  the  structural  stiffness  matrix  K  of  Eq .  (4) 
(2).   The  effects  of  shear  deformation  may  be  discarded  by 
setting  r  =  0  in  Eq .  (8). 

3.2.   Certain  Effects  Related  to  Pressure,  Rigid  Parts  and 
Supports.   The  generalized  coordinates  are  taken  to  be  the  dis- 
placements y.,  the  rotations  9.,  the  displacement  E,    of  the 
projectile  with  respect  to  the  breech,  and  the  axial  displacement 
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u  of  the  barrel  and  the  breech  (as  noted  previously;  see  Fig.  1). 

The  breech  is  regarded  as  a  rigid  block  of  mass  M  and  moment  of 

inertia  I  about  its  center  of  mass.   The  center  of  mass  of  the 

breech  may  be  off  center  from  the  axis  of  the  barrel.   The 

eccentricity  of  the  breech  is  specified  by  two  lengths,  d  and  e 

(Fig.  1).   There  is  an  axial  spring  with  constant  aR  and  an 

axial  dashpot  with  constant  u_  at  the  breech.   These  elements 

are  offset  distances  er  and  e"  from  the  axis  of  the  barrel  (Fig.  1) 

Movement  of  the  breech  transverse  to  the  barrel  is  restrained 

by  a  spring  with  constant  a-  and  a  dashpot  with  constant  u, . 

These  are  included  among  springs  and  dashpots  with  constants 

a.  and  u.  that  act  at  the  nodal  points  on  the  barrel.   Their 
J       J 

effects  will  be  discussed  later. 

3.2.1.   Action  of  Gravity  on  the  Breech.   By  Fig.  1,  the 
descent  of  the  center  of  mass  of  the  breech  due  to  displacement 
at  x  =  0  is 

u  sina  -  y,  cosa  +(d  cosa  -  e  sina)9, 

where  a  is  the  angle  of  inclination  of  the  barrel  with  respect 
to  the  horizontal.   Therefore,  the  contribution  of  the  action 
of  gravity  on  the  breech  to  the  variation  of  the  virtual  work  is 

Mglfiu  sina  -  5y,cosa  +(.d  cosa  -  e  sina)59,]  (9) 

By  Eqs .  (2}    and  (3),  it  is  seen  that  the  coefficients  of  5u, 
<Sy-.,  <59,  in  Eq.  C9 )  contribute  to  the  right-hand  side  of  Eqs.  (3). 
With  Eqs.  C3 )  written  in  the  form  of  Eqs.  (4),  the  coefficients 
of  5u,  <5y,,  59-.  contribute  to  the  right-hand  side  of  Eqs.  (4). 
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3.2.2.  Contribution  to  6W  of  Axial  Force  on  Projectile. 

2 
The  projectile  is  subjected  to  the  axial  force  ira  p-F,  where  p 

is  the  gas  pressure  behind  the  projectile,  a  is  the  radius  of 

the  bore,  and  F  is  the  axial  frictional  force  of  the  barrel  on 

the  projectile.   The  effect  of  the  shock  wave  ahead  of  the 

projectile  can  be  included  in  F.   The  contribution  to  5W  of  the 

axial  force  on  the  projectile  is,  therefore, 

(^a2p-F)(5£  -  6u)  (10) 

3.2.3.  Contribution  to  <5W  of  Axial  Movement  of  the  Barrel. 
Since  pressure  p  also  acts  on  the  breech,  the  contribution  of 
the  axial  movement  of  the  barrel  to  6W  is 

2 

Mg(sina)5u  +  Oa  p-F)6u  (11) 

where  M  is  the  mass  of  the  barrel,  and  where  the  term  Mg(sina) 
is  the  contribution  of  the  axial  component  of  gravity. 

Accordingly,  the  net  contribution  to  5W  from  the  breech, 
the  axial  force  acting  on  the  projectile,  the  axial  movement 
of  the  barrel,  the  breech  spring,  and  the  breech  damper  is, 
by  Eqs.  (9),  (10),  (11), 

6¥  =  [-aBCu-ef81)-uB(u-e"9]_) 


(12) 


+  CM+M)g  sina]5u  +  (7Ta2p-F)5? 
-Mg(cosa)6y,  +  [Mg(d  cosa  -  e  sina) 
+aBe'(u-e'81)+  uBe" (u-e"9 ± ) J  5 S± 

The  terms  involving  aB  and  y0  in  Eq.  (15)  come  from  the  fact 

that  the  compression  of  the  breech  spring  is  u-ef9,.   Consequently, 


the  resisting  forces  of  these  elements  are  a  (u-e'0,)  and 
U_(u-e"9,).   The  coefficients  of  <5u,  6£  ,  5y,  ,  60,  contribute 
to  the  corresponding  terms  in  Eqs .  (3)  or  (4). 

3.2.4.   Contribution  to  6W  of  Dampers  and  Supporting 
Structure.   The  contribution  to  6W  from  rotational  and  transla- 
tional  dashpots  at  point  x  .  is 

-y  .y  .6y  .  -  u'.9  -50  .  (13) 

J  J  3  3    3      3 

where  ( y  .  ,  u '.  )  are  constants  for  the  dashpots,  and  the  dot 

J    J 

denotes  time  derivative.   At  a  node  where  there  is  no  dashpot , 

y  .  =  u  '.  =  0  . 
J     J 

If  the  supporting  structure  can  be  regarded  as  a  set  of 
decoupled  springs  located  at  the  nodes,  the  contribution  of  the 
springs  at  node  j  to  6W  is 

-a.y  .<Sy  .  -  a*.  9  .59  .  (14) 

J  J   J    J  J   J 

Accordingly,  the  contributions  of  the  supports  and  dampers  to 
6¥  is 


n+1  n+1 

6W  =  -  Y^  (u  .y  .+a  .y  .  )6y  .  -  V^  (y  \B  .+a  r.9  .  )69  .         (15) 

J-i  j-i 


Again,  the  coefficients  of  6y  .  and  50.  contribute  in  the  appro- 
priate  manner  to  Eqs.  (3)  and  (4).   In  particular,  the  coeffi- 
cients of  y . ,  0 .  contribute  to  the  matrix  K  and  the  coefficients 
of  y,  0.  contribute  to  the  matrix  C  in  Eq .  (4). 


19 


3.2.5.   Action  of  Gravity  on  the  Barrel.   In  the  finite 
element  method  employed,  the  external  loads  (forces  and  couples) 
are  applied  at  the  nodes  only.   If  a  distributed  load  or  a  con- 
centrated load  is  applied  between  nodes,  they  must  be  replaced 
by  equivalent  systems  of  concentrated  loads  at  the  nodes.   We 
require  these  equivalent  loads  to  be  derived  in  such  a  manner 
as  to  be  consistent  with  the  method  used  in  deriving  the  stiff- 
ness matrix  (Appendix  A).   Since  the  lateral  component  of  the 
gravity  load  acting  on  an  element  of  the  barrel  is  uniformly 
distributed  between  nodes,  the  equivalent  nodal  force  matrix 
for  the  jth  element  of  the  barrel  is 


-m.A.  g  cos a 

a    =      J    3 

i         12 


,  j  =  1,2, n   (16) 


where  m.  =  the  mass  of  the  jth  element  and  A.  =  the  length  of 
J  J 

the  jth  element.   By  the  method  of  finite  elements,  appropriate 
additions  of  the  matrices  f  .   are  required  to  incorporate  the 
element  nodal  forces  into  the  driving  force  vector  f(t),  Eq .  (4) 

The  effect  of  the  axial  component  of  gravity  on  the  barrel 
has  been  accounted  for  previously  (see  Eq .  11). 

3.2.6.   Action  of  Gravity  on  the  Tuning  Mass.   In  general, 
the  tuning  mass  is  located  at  some  point  between  nodes,  Figs.  1 
and  2.   Since  the  lateral  component  of  gravity  force  on  the 
tuning  mass  is  mT  g  cosa ,  this  force  must  be  transformed  into 
equivalent  nodal  forces  for  the  element  in  which  the  tuning  mass 
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is  located.   By  the  method  of  finite  elements  (Appendix  A),  the 
equivalent  nodal  force  matrix  for  the  element  in  which  the  tuning 
element  lies  is 


-  3s2  +  2s3 


fT  =  -mT  g  cosa 


(1    -  2s   +   s  )sA 


3s' 


2s' 


(_-s   +   s  )sA 


C17) 


where  A~  =  the  length  of  the  tuning  mass  element  and  sTAT  =  the 
distance  from  the  left  node  of  the  element  to  the  tuning  mass. 
Equation  (17)  accounts  for  the  lateral  component  of  gravity 
acting  on  the  tuning  mass.   In  addition,  the  axial  component  of 
gravity  contributes  to  <5W  in  the  form 


5¥  =  mT  g  (sina)6u 


(18) 


3.2.7.   Action  of  Gravity  on  the  Projectile.   The  absolute 
axial  displacement  of  the  projectile  at  time  t  is  (5-u),  Figs. 
1  and  3.   Then,  the  contribution  to  5W  of  the  axial  component 
of  gravity  on  the  projectile  is 


6W  =  mg  (_sina)  (5u  -  6?) 


(19) 


The  effect  of  the  lateral  component  of  gravity  on  the  projectile 
is  included  in  the  same  manner  as  for  the  tuning  mass  (Eq.  17) 
where  A? ,  sT  are  replaced  by  ApJ  s    respectively.   In  addition, 
it  is  to  be  noted  that  s   is  a  function  of  time  t,  since  %    =  £(t) 
In  studying  the  effect  of  gravity  on  the  projectile,  it  has  been 
found  by  actual  computation  that  this  effect  is  extremely  small 
for  projectiles  of  the  20-40  mm  class. 
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3.2.8.   Axial  Friction  Effect  of  Projectile  on  Barrel.   The 
axial  frictional  force  F(t)  exerted  by  the  projectile  on  the 
barrel  is  assumed  to  be  a  known  function  of  time  t.   Alternatively, 
it  may  be  computed  if  the  pressure  acting  on  the  projectile  and 
the  projectile  acceleration  are  known  as  functions  of  time. 
Also,  it  may  be  expressed  as  F(£),  where  x  =  C(t)  locates  the 
projectile  in  the  barrel.   This  is  a  distance  s  along  the  curve 
of  the  axis  of  the  barrel  (Fig.  4). 

There  is  a  subtlety  in  the  calculation  of  the  work  that 
force  F  performs  on  the  barrel,  because  the  projectile  is 
sliding  in  the  barrel.   For  example,  if  a  grinding  wheel  acts 
on  a  fixed  plate,  the  frictional  force  of  the  grinder  performs 
no  work  on  the  plate  because  the  plate  does  not  move.   However, 
the  frictional  force  of  the  plate  performs  negative  work  on  the 
grinder. 

The  axial  movement  of  the  projectile  during  an  infinitesimal 
time  dt  does  not  affect  the  work  that  the  projectile  performs 
on  the  barrel.   Consequently,  in  computing  5W,  we  set  dt  =  0. 
The  corresponding  point  on  the  curve  y  +  5y  consequently  also 
lies  a  distance  s  along  the  curve.   The  component  of  F  along  the 
virtual  displacement  determines  5W  (Fig.  4).   We  have 


Xl 


y  =  y(x),  s  =  /     /l  +  (y» )z    dx 

o 


Also,  along  the  curve  y  +  5y, 


X2 


s  =  /     /l  +  ( y  '  +  5  y  '  ) 2  dx 
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Since  y  <<  1 , 


r1 

s  «  x±   +  I  /   (y')2  cLx 


and 


r2 

s  «  x2  +  I  /   (y'2  +  2y'6y>)  dx 


Therefore, 


Xl  *2 

x2"xi  =  £  =  \  I     (y')2dx  -  |  /     (y'2  +  2y,<sy')  dx 


and  since  x,  =  E,  ( t )  , 


£  =  /   y'2  dx  -  |  /    (y'2  +  2y'6y'  )  dx 
o  o 


and  we  have 


£  -  \    I     y,2dx-^/  (y'2  +  2y'5y' )  dx 
o  o 

-  /   y'2  dx 
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or  r  f, 

S  ?  +  £ 


=  -/  y'Sy'dx  -  |  /    y'2  dx 


5 
However,  since  £  is  small, 


/ 


y'2  dx  =  £[y'(?)]2 


Hence , 


£{i  +  |[y'(?)]2}  =  -/  y'Sy'dx 


and  since  y'  <<  1 

? 

£  a  -  /   y  '  5  y  '  dx 
o 

Using  the  trapezoid  rule  to  evaluate  this  integral,  and  setting 
E,   =  x  ,  we  get 

e  =  -  i[(8169lte2Je2)A1+(e2se2+936e3)A2 
+  (936e3+946e4)i3  +  .... 

+  (6r-l59r-l  +WAr-l] 

The  contribution  to  5W  from  axial  friction  is  the  component 
of  F  in  the  x  direction  times  £  plus  the  component  of  F  in  the 
y  direction  times  Sy  .   Thus,  we  have  approximately 

5W  =  F£  +  F9  Sy 
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Accordingly,  for  the  j-th  interval, 

6W.  =  -fo^+e^se^Aj,  i  <  3   <  r-i 

5W.  =  0  ,  j  >  r 
J         ~ 

In  addition,  there  is  the  contribution  F0  6y   at  point  x  . 

Consequently,  the  contribution  of  axial  friction  of  the 
projectile  on  the  barrel  to  6W  is 

6W  =  F9r0yr  -  \   FA191661  -  \   F(  A-L+A2  )92692 


(20) 


"  \   F(A2+A3)93693  " 


-  =?   F(A   0+A   ,  )8   ,59   , 
2   s  r-2   r-l;  r-1   r-1 

-  I  F  A   ,9  69   , 
2     r-1  r   r  ' 


As  with  previous  contributions  to  6W,  the  terms  in  Eq .  (20)  must 
be  incorporated  appropriately  into  Eqs.  (3)  or  (4),  with  the 
understanding  that  F  =  F(t),  or  alternatively,  F  =  F(£).   In 
particular,  terms  in  Eq.  (20)  contribute  to  the  K  matrix  of 
Eq.  (4). 

3.2.9.   The  Bourdon  Effect.   The  tendency  of  the  pressure 
p  in  a  gun  barrel  to  straighten  the  barrel  is  called  the 
"Bourdon  effect".   A  line  element  of  a  generator  of  the  bore 
above  the  neutral  plane  has  length  (1  -  ay    sin9)dx  (Fig.  5). 

XX 

The  corresponding  element  below  the  neutral  plane  has  length 

(1  +  ay   sin9)dx.   The  difference  of  these  lengths  is 

2ay   sin9  dx ,  and  the  difference  between  the  areas  of  infini- 

XX 

2 

tesimal  strips  is  accordingly  2a  y    sin9  d0  dx.   The  projection 

XX 

9        2 
of  this  area  on  the  neutral  plane  is  2a  y    sin  9  d9  dx . 

XX 
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Consequently,  the  effective  load  per  unit  length  of  the  barrel, 
caused  by  the  gas  pressure  is 


IT 

2/2  2 

-  2a  py    |    sin  9  d0  =  -  ira  py 

XX  I  XX 

o 


(21) 


Accordingly,  the  contribution  of  the  gas  pressure  to  5W  is 


6W  =  -  Tra  p 


y    6y  dx 
xx  J 


0 


where  E,    is  the  location  of  the  projectile.   We  suppose  that  the 
projectile  lies  at  a  net  point  x  ;  i.e.  ,  E,   =  x  .   This  is  not 
unrealistic,  since  the  terminal  point  of  the  pressure  is  rather 
indefinite  because  of  the  shock  wave  ahead  of  the  projectile 
and  the  finite  length  of  the  projectile.   By  Eq.  (21), 


x. 


5W  =  -  na  p 


y    oydx  +   /    y    5ydx  + 

J XX   J         /       XX   J 


(22) 


Hence , 


x 


r-1 


yxx  6ydx 


r-l 

6W  =  2_j      6W  • 

J  =  l 
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where  6W .  is  the  contribution  of  the  "Bourdon  effect"  to  5W  in 

the  j-th  interval  (A   =  x   ,  -  x  ).  We  let  <5W   =  0  if  j  >  r-1 , 
since,  for  x  >  x  ,  p  =  0. 

In  the  j-th  interval,  Eq .  (1)  of  Appendix  B  yields 

y        =    2a^    +   6a^    x 
J  xx  2  3 

and  Sy   =    6aj5)   +   x6a^    +   x      6a^    +   x      6a^ 


Hence , 


Xj+1 


5W.    =    -TTa2p    /  (2a^    +   6a^x)(6ajj   +   x6a^    +   x25a^ 


x  . 

3 


+   x36a^)    dx 


Integration  yields 

5W      =    -Tra2p{2(xj+1    -    x.)    a^5aj 


+    (x2  +  1    -   x2)(aJ5aJ    +   3a^a£) 


(23) 


+  |(xo5+1    -   Xj5)    a^} 

Introducing  Eq.  (2)  of  Appendix  B  into  Eq .  (23),  we  get  an 
expression  of  the  following  form: 
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6wj  -  BiiVyj +  Bi2^+i^j +  Bi36j^j +  Bi49J+1^j 

+   B21yj6yj+1    +   B22yj+l5yj+l    +   B23Vyj+l    +   B249j+l^j  +  l 
+   B3lV9J    +   B3VJ  +  lSej    +   BJ33ejSej    +   B349j  +  lMj         <24> 
+   B4lWl    +   B42yj+lSVl    +   B43Wl    +   B449J+l"j+l 

If   we    introduce   Eq .     (2)    of    Appendix   B    into   Eq .     (23),    and    compare 
the    resulting   lengthy    expression   with   Eq .     (24),    we    get    expressions 
for    the    coefficients    BJ _ .       With   Eq .     (3)    of    Appendix   B,    they 
simplify    greatly.       Finally,    they    reduce   to    the    following: 

2  2 

rJ       =    r!       =    6jtsl    P  rJ       =    rJ       =         6ira"p 

Dll         D22            5A  .  '                         12         D21                 5A  . 

J  3 

2  2 

j      _    llTra    p  j  Tta   p 

D13  10         '  D31  10 

2  2  2 

oJ  =     oJ  =     ^a     P  Rj  =  UTTa     P  Rj  TTa     p 

14         D41            10       '                         D24                   10  '    D42  10 

2                             .                          27Ta2pA  . 
B23    "    B32 10"    '  B33    "    B44 1~5  (25) 


2 

M 
'34         "43  30 


BJ       =    BJ       =    -   ^-^ 


These  equations  hold  for  j  >  r-1,  with  p  =  0.   Hence,  B^ „  =  0 

if  j  >  r-1.   Also,  the  equations  hold  for  r  =  1;  j  =  r-1  =  0  with 

the  understanding  that  B  n    =  0.   The  virtual  displacements  5u 

to       aft 

and  SE,    contribute  nothing  to  the  Bourdon  effect.   It  is  to  be 

noted  that  the  coefficients  B  n    are  not  entirely  symmetrical  in 

ag 

a  and  B.   The  terms  in  Eq .  (24)  contribute  in  an  appropriate  manner 
to  the  stiffness  matrix  of  Eq.  (4). 
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3.2.10.   Bending  Caused  by  Axial  Inertia.   The  tension 
P(x,  t)  in  the  barrel  due  to  recoil  acceleration  may  be  calcu- 
lated by  xNTewton's  law  with  the  assumption  that  the  barrel  is 
rigid.   Accordingly, 

P  =  mil  (26) 

where  m    is  the  mass  of  the  part  of  the  barrel  in  the  range 
greater  than  x. 

If  we  suppose  that  the  barrel  is  inextensional ,  the  shorten- 
ing of  the  chord  due  to  bending  is 


I  /   yx  dx 


where  L  =  x  ,,  is  the  length  of  the  barrel.   Effectively,  then, 
n+1 

1   2 
an  element  of  length  dx  is  shortened  by  the  amount  ~-  y  dx ,  and 

1    2 
the  potential  energy  of  that  element  due  to  shortening  is  -~   Py  dx 

Accordingly,  the  strain  energy  is  augmented  by  the  amount 

L 

\  f      P?x  <* 
0 


The  integral  is  a  function  of  t,  since  y  and  P  are  functions 
of  x  and  t.   Section  properties  of  the  barrel  do  not  enter  here, 
except  insofar  as  they  affect  P(x,  t\.      Equation  (1)  of  Appendix  B 
now  yields 

x  . 
V  .  =  j        I  PCa^  +  2a^x+  3a^  x2)2  dx        (.27 )_ 
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in  which  V.  is  the  potential  energy  in  the  j-th  interval  due  to 
shortening  of  the  chord. 

We  may  write  V   in  the  following  form,  by  virtue  of  Eq .  (2) 
of  Appendix  B, 

2V.    -   AJlVj    ♦  AJ2yj  +  1yj+   A^e^    ♦   A^B^y. 

+   A21yjyj+1    +   A22yj+lyj+l    +   A238jyj+1    +   A249j+lyj+l 

♦    4^9.    ♦    AJ2yj  +  l9j    +    A^ejBj    +    AJ49j+l9j  (28) 

+  Aiiyj9J+i +  Ai2yj+iaJ+i +  Ai39j9J+i +  Ai4ViVi 

Without  loss  of  generality,  we  may  set  AJ _  =  A~  . 

It  is  reasonable  to  assume  that  P  varies  1 inearly  in  an 

interval.   Then,  in  the  ith  interval  x.  <  x  <  x.^-,) 

3  J  +  l' 

P  =  P.A"1(x.J.1  -  x)  +  P.^,A-1(x  -  x.)  (29) 

3    3         J  +  l  J  +  l  J        3 

Consequently , 


H+1  P 

/     n,,.  i        r  n+2  ,  lON  n  +  1     .,  .,  .  n  +  2-, 

/    *  Pdx  =  (n+l)(n+2)aj  txJ  +  1-(^2)Xj   xj  +  1  +  (n+l)Xj 

x.  (30) 

3  p 

i+1  r  n+2  ,   „,    n+1   ,   ,  .  n+2, 

+  (n+l)(n,2)Aj  {xj   -(n+2)xjxj+l  +tn+l)xJ+1} 

With  Eq.  (29),  the  integrals  in  Eq .  (27)  can  be  evaluated.   Then 

with  Eq.  (2)  of  Appendix  B,  we  get  V.  expressed  as  a  quadratic 

form  in  y . ,  y..n,  9.,  8..,.   Using  Eq .  (3)  of  Appendix  B,  and 
J    J+l    3         J+l 

comparing  with  Eq .  (28),  we  get 
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All  ■  A22  "  I  AJ1(Pj+Pj+l) 

A33  =  55  &J1(3Pj+PJ+lJ"  A44  -  JO   tfl*}**!^              (31) 

A12  =  A21  =  "  5  Aj  ^j^j+l*'   A13  =  A31  =  TO  Pj  +  1 

aJ   =  a-3   =  1  p  aJ  =  a^  -  -  — -  P 
14     41    10  y                                   23     32      10   j+1 

A24  "  A42  "  "  TO"  Pj'  A34  "  A43  "  "  W   W'j  +  l' 

A°K  =  0,  An*l=  0 


The  contribution  to  5W  due  to  bending  caused  by  axial 
inertia  is 

6Wj  =  -A^y^yj  -  AJ2y.  +  16yj  -  A^e^  -  A^^ 

~4l*  1**3+1    -  *aaWj+l  "  A239jSyJ+l  "  A249j+l6yj+l 
-A31yjS9j  "  *32yJ+l"j  "  A33ej59j  -  A349J+l69j       <32> 
-A4iyj69j+1  "  A42yj+l69j+l  "  A439j59j+1  "  A446j+l69j+l 

The  terms  in  Eq.  (32)  contribute  in  an  appropriate  manner  to 
Eq.  (4). 

4.   KINETIC  ENERGY  EFFECTS 

The  kinetic  energy  is  the  sum  of  the  kinetic  energies  of 
the  breech,  the  tuning  mass,  the  projectile,  and  the  barrel. 
These  energies  contribute  to  the  mass  matrix  of  Eq.  (4)  or  the 
total  kinetic  energy  T  of  Eq .  (3). 
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4.1.   Kinetic  Energy  Contribution  of  the  Breech.   The  veloc- 
ity components  of  the  center  of  mass  of  the  breech  are  e9,  -  u, 
y,  -  d6,  (see  Fig.  1)  and  the  angular  velocity  of  the  breech 
is  9,.   Consequently,  the  kinetic  energy  of  the  breech  is 

Tb  =  1  M(e61  -  u)2  +  \   M(yi  -  d^)2  +  |  I^2 


where  M  is  the  mass  of  the  breech  and  I  is  the  moment  of  inertia 
of  the  breech  about  a  transverse  axis  through  its  center  of  mass. 

Accordingly  by  Eq .  (3),  the  breech  contributes  the  term 

-  ..  --22 

M(y   -  9  d)  to  the  first  row  of  Mx  and  the  term  [I+M(e  +d  )]• 

9,  -  Mdy,  to  the  second  row  of  Mx,  and  in  addition,  the  term 

Meii  must  be  added  to  the  right  hand  side  of  Eq .  (4)  in  the  second 

row. 

4.2.   Kinetic  Energy  Contribution  of  the  Tuning  Mass.   The 
contribution  of  the  tuning  mass  to  the  kinetic  energy  may  be 
treated  by  means  of  D'Alamberts  Principle.   Thus,  we  represent 
the  kinetic  effects  of  the  tuning  mass  by  an  inertial  force  m„y 
and  an  inertial  couple  I_9  directed  as  shown  in  Fig.  2.   Hence 
representation  of  y  (9  =  y')  by  a  cubic  for  the  element  in  which 
the  tuning  mass  lies  (see  Eq .  1,  Appendix  A  or  Eq.  1,  Appendix  B) 
yields  the  concentrated  force  mTy  and  the  concentrated  couple 
I„9  as  functions   of  the  nodal  displacements  and  nodal  rotations 
of  the  tuning  mass  element  ,  x_   x  <_   x„,  ,  .   Then  by  the  results 
for  a  concentrated  force  and  a  concentrated  couple  (see  Appen- 
dix A)  mTy  and  IT9  may  be  transformed  into  equivalent  nodal 
forces  and  couples. 
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4.3.  Kinetic  Energy  Contribution  of  the  Barrel.   The  kinetic 
energy  of  the  barrel  contributes  to  the  term  Mx  in  Eq.  (4).   By 
the  method  of  finite  elements,  the  mass  matrix  m    for  a  beam 
element  may  be  computed  and  the  contributions  to  the  total  mass 
matrix  M  may  be  obtained.   (See  Appendix  A  for  the  mass  matrix 
m.) 

4.4.  Kinetic  Energy  Contribution  of  the  Projectile.   As  for 
the  barrel  element,  the  projectile  kinetic  energy  contributes 

to  the  term  Mx.   However,  in  addition  the  projectile  kinetic 
energy  contributes  to  the  terms  Cx  and  Kx ,  because  of  the  nature 
of  its  motion  (Coriolis  acceleration,  centripetal  acceleration, 
etc.).   The  analysis  of  these  contributions  are  best  undertaken 
by  a  direct  consideration  of  the  acceleration  effects  of  the 
projectile  upon  the  equations  of  motion  of  the  system  (see 
Appendix  C) . 

5.   THE  COMPUTER  PROGRAM  AND  RESULTS 

The  computer  program  uses  a  finite  element  method  to  analyze 
the  dynamic  motion  of  a  gun  barrel.   The  problem  is  illustrated 
in  Fig.  6.   The  projectile  tuning  mass  and  breech  are  considered 
to  be  rigid  bodies.   The  barrel  is  treated  as  a  beam  and  is 
divided  into  a  number  of  finite  elements.   The  division  points 
for  these  elements  are  called  nodes  or  node  points,  and  there 
are  two  degrees  of  freedom  at  each  node;  lateral  deflection  y, 
positive  upward  and  rotation  y',  positive  counter  clockwise. 
The  nodes  are  numbered  from  the  left  end  (the  breech  end)  towards 
the  muzzle  (see  Fig.  7).   If  there  are  N  elements,  there  are  N+l 
nodes  and  2N+2  degrees  of  freedom. 


33 


Each  of  the  elements  may  have  a  different  length.   The  barrel 
may  be  subdivided  into  a  number  of  segments.   Within  each  segment, 
the  elements  are  the  same  length,  but  they  may  differ  in  diameter 
if  the  barrel  is  tapered  (Fig.  7).   The  tapered  portions  are 
modeled  in  a  "stairstep"  manner. 

The  individual  elements  are  treated  as  beams  with  four 
degrees  of  freedom  (deflection  and  rotation  at  both  ends). 
Usually,  the  Bernoulli-Euler  beam  theory  is  sufficient  for  most 
problems.   However,  the  program  provides  the  option  of  including 
shearing  deformation  and  rotary  inertia  (Timoshenko  beam  theory). 
In  order  to  include  shearing  deformation,  a  cross-section  shear 
constant  must  be  input.   This  may  be  obtained  from  the  literature 
[3].   For  a  hollow  tube  cross  section,  the  value  of  the  shear 
constant  is 

SHEAR  =  6(l+v)(l+m2)2  _ 


(7+6v)(l+m2)2  +  (20+12v)  m2 


,     -    inside  radius      „  .     , 

where  m  =  — .■.  , 3-r — ,  v  =  Poisson  s  ratio 

outside  radius' 


The  program  permits  translat ional  and  rotational  springs 
and  dashpot  supports  at  any  number  of  node  points  in  addition 
to  a  recoil  spring  and  dashpot  at  the  breech.   A  Bourdon  pressure 
and  a  projectile-barrel  friction  force  are  incorporated  in  the 
program,  and  may  be  included  at  the  option  of  the  user. 

The  program  is  intended  primarily  to  be  used  for  a  numerical 
time  integration  study  of  a  gun  barrel  motion  from  the  time  of 
first  firing  onward  including  the  option  of  more  than  one  round. 
The  numerical  integration  routine  utilizes  the  Newmark  "Beta 
Method"  [4].   If  desired,  the  program  can  be  used  for  determination 
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of  the  static  deflected  shape  with  no  time  variance  involved. 
In  addition,  provision  for  a  shaking  of  the  supports  of  the 
barrel  (due  to  vehicle  motion  or  other  cause)  is  available  as 
an  option. 

The  program  output  is  printed  as  follows:   (1)  the  problem 
title,  (2)  the  input  data,  (3)  the  initial  values  of  the  posi- 
tion, velocity  and  acceleration  of  the  node  points,  (4)  the 
position,  velocity,  and  acceleration  of  the  node  points  at  time 
intervals  dictated  by  the  input  requests,  (5)  a  plot  of  the 
muzzle  deflection  as  a  function  of  time. 

Typical  plots  of  muzzle  deflection  as  a  function  of  time 
are  shown  in  Figs.  8-16.   The  data  is  for  the  British  Rarden 
system  for  a  range  of  spring  stiffnesses  K,  damping  coefficients 
C  and  various  locations  of  barrel  lateral  stopper  pads.   Figure  8 
is  for  a  single  shot,  whereas  Figs.  9-16  are  for  bursts  of  three 
shots.   Figures  8-9  are  for  shots  without  the  lateral  motion 
stopper  pads,  and  Figs.  10-16  are  with  pads  at  various  locations 
of  the  pads.   Study  of  these  figures  reveals  the  great  importance 
of  the  stopper  pads.   In  general,  the  use  of  a  stopper  pad 
reduces  the  tip  (muzzle)  deflections  by  a  factor  of  ten  or  more. 
Study  of  Figs.  11-14,  for  the  case  of  the  pad  located  54  inches 
from  the  breech  end  of  the  barrel,  shows  that  increasing  the 
damping  coefficient  C  also  tends  to  reduce  the  muzzle  deflection, 
particularly  in  the  later  firings.   However,  a  predominant  role 
is  played  by  the  stopper  pad  in  the  control  of  muzzle  displacement 

With  appropriate  input  data,  other  gun  systems  may  be  readily 
studied , 
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APPENDIX  A:   FINITE  ELEMENT  MODELS 

Al .   Beam  Element 

A  typical  element  of  beam  has  4  degrees  of  freedom  as  shown 
in  Figure  Al .  The  displacement  of  a  typical  point  in  the 
element  is  assumed  to  be 


y  =  aQ  +  a,x  +  a2x   +  a„x 


(1) 


W 


ith  the  four  boundary  conditions  y  (  0  )=u-,  ,  y  '  C01=Uo  ,  y(_A  1=^3  »  y  f  CA1= 


u.  ,  the  displacement  at  a  typical  point  can  be  expressed  in  terms 
of  the  nodal  displacements  (u, ,  u« ,  u„,  u.)  as  follows: 


r  -,  _  r ,    3x     2x       2x     x    3x 

[y]  "  [     "72    —3'  x"  ~  +  "2'  ~2 
A      A  A     A 


or 


y  =  Nu 


o  3      2     3 
2x      x     x  -I 

A3'  "  A     A2 


1 —            — 

ul 

U2 

U3 

U4 



(2) 
(2a) 


The  slope  y'  may  thus  be  expressed  in  terms  of  the  nodal  dis- 
placement as 


r  ,  -i    r_  ^2i  +  6x    i    4x    3x    6x    6x      2x    3x 
LY         A2     A3'       *     A2'  A2  "   A3'  ~~         "? 


—            — 

ul 

U2 

U3 

U4 

1 

(3) 
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The  curvature  y"  may  likewise  be  expressed  in  terms  of  the 
nodal  displacements  as 


[y..]  =  [-  «   +  i2| 

A^     AJ 


4    6x   6_    12x     2    6x 

"  A    A2'  A2  "   A3'  "  A    A2' 


— 

— 

ul 

U2 

U3 

U4 





(4) 


£  =  Bu 

The  relation  between  moment  and  curvature  is  (ignoring  shear 
deformation;  see  Art.  A2 . ) 


[M]  =  [EI]  [y"] 


(5) 


a  =  Ee 

The  general  expression  for  stiffness  matrix  is  (derived  by  the 
principle  of  virtual  work)(2) 


k  = 


B   EB  dx  =  EI 


12 

3 


A 


A 
12 

A3 

_6 
,2 


12 

A3 


A 

12 
,3 


A 


(6) 
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This  stiffness  matrix  for  the  beam  element  relates  nodal  loads 
to  nodal  displacements.   The  nodal  loads  are  shown  in  Fig.  A2 . 


'11     *12 


—     — 

ul 

u2 

U3 

U4 



(7) 


f  =  ku 

The  equivalent  load  matrix  for  a  distributed  load  is 
obtained  from 


eq 


p  dx 


(8) 


where  p  =  [p]  is  the  distributed  load,  e.g. ,  if  p  =  p   (a  con- 
stant )  ,  then 


:eq 


N   [p  ]  dx  =  p 
^o        *o 


1/2 
A/12 
1/2 
-A/12 

(1)   Uniform  distribution  load  equivalent 


(9) 


wA 
2 


w 


>;    t v   ±      \> 


wA 
12 
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wA 
2 

) 

A2 
wA 

12 


(2)   Linear  distributed  load  equivalent 


c 


3wA 
20 


wA 
30 


7wA 


20 


5 


wA'" 
"20" 


(3)   Quadratic  distributed  load  equivalent 


{ 


wA 
15 


wA 
"60" 


4wA 


15 


) 


wA' 
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If  a  concentrated  load  is  to  be  applied  at  some  point  other 
than  a  node,  it  must  be  replaced  by  an  equivalent  system  of  nodal 
loads  also.   Examples  of  equivalent  systems  are  given  below: 
(4 )   Concentrated  force 


F 
sA  -» 


F(l-3s2+2s3) 


FsA(l-2s+s  ) 


2    3 
F(3s  -2s  ) 


FsA(s-s") 
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(5)   Concentrated  couple 


K  sA  -| 


Ci7 


A 


6c  r        2,  6c  ,        2., 


A 


( 


A 


A 


2 


cCl-4s+3s2l  c(;2sr-3s2) 


A2 .   Beam  Element  Including  Shear  Deformation 

Consider  a  beam  element  of  length  AZ  (Fig.  A3).   Under 
deformation  cross  section  plane   TOB  of  the  beam,  which  is 
initially  perpendicular  to  the  beam  axis,  rotates  into  the 
position  T*0*B*.   Because  of  shear  effects,  the  plane  T*0*B* 
is  no  longer  perpendicular  to  the  displaced  axis  of  the  beam, 
Fig.  A3,  but  rotates  through  an  angle  tj>  f    y'.   For  equilibrium 
of  an  infinitesimal  length  dZ  of  the  beam  element,  Fig.  A4 , 
we  have 


dM        dF    _ 
dZ  '     '   dZ     p 


(10) 


where  M  denotes  the  bending  moment ,  F  denotes  the  shear  force 
and  p  denotes  the  lateral  pressure  (+  as  shown)  that  act  on  the 
element  dZ .   The  stress-strain  relations  for  the  beam  element 
are 


M  =  EI<J>'  ,   F  =  GA(y  • -<J> ) 


(11) 


where  E  =  the  modulus  of  elasticity,  G  =  the  shear  modulus,  I  = 
the  moment  of  inertia  of  the  area,  A  =  the  beam  cross  sectional 
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area,  and  prime  denotes  derivative  with  respect  to  axial  length 
Z. 

Substitution  of  Eq .  (11)  into  Eq .  (10)  yields  the  equilibrium 
equations  in  terms  of  $    and  y  as 


EI<|)M  =  GA(<J)-y  *  ) 
GA(yM-<J>»  )  =  -p 


(12) 


In  the  simplest  finite  element  shear  model,  the  displacements 
are  assumed  in  the  form  (2) 


cj>  =  A1   +  A2Z  +  A3Z 

y  =  A4  +  A5Z  +  AgZ2  +  A?Z3 


(13) 


Assuming  that  the  pressure  p  (or  any  other  load  between  node 

points  of  the  beam  element)  is  transformed  into  nodal  forces 

(see  Eq.  3,  Art.  Al ) ,  we  find  that  substitution  of  Eqs .  (13) 
into  Eqs.  (12)  yields  the  result 


<$>    =  Ax  +  A2Z  +  A3Z' 


y  -  A4  +  CA,  -^  VZ  +  ^z2  +^z3 


(14) 


where 


r  = 


12EI 
AGA2 


(15) 


is  the  shear  deformation  ratio.   With  Eqs.  (14)  the  method  of 
finite  elements  (2)  yields  the  beam  element  stiffness  matrix 


k  = 


EI 


A  (1+r) 


12    6A 


-12 


6A 


6A    A2(4+r)  -6A     A2(2-r) 


-12   -6A 


12    -6A 


6A    A2(2-r)  -6A     A2(4+r) 


(16) 
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We  observe  that  if  r  =  0,  Eq.  (16)  reduces  to  Eq.  (6). 
The  mass  matrix  m    for  the  beam  element  is  given  by 


m   = 


N   p  N  dZ 


(17) 


where  the  N  matrix  relates  the  end  (nodal)  displacements  x, ,  x„ , 
x„ ,  x.  of  the  beam  element  to  the  displacements  y,  <j>  (see  Eq .  2) 
and  the  p  matrix  is  the  mass  density  matrix  (2) 


P  = 


UR 


(18) 


where  u  =  mass  per  unit  length  and  R  =  the  radius  of  gyration 
of  the  beam  cross  section. 

Explicitly,  the  4x4  matrix  m   may  be  written  as 


m    = 


PA 
840(l+r): 


m 


.77 


m 


m 


11 


12 


13 


14 


71 


12 


'77 


777 


,77 


22 


23 


24 


-77 


,77 


.77 


77 


13 


23 


33 


34 


m 


77 


14 


24 


77 


34 


,77 


(19) 


44 


43 


where 

2 

w.n1  =  312  +  588r  +  280r2  +  1008  ^5- 

2 

m12   =  A[44  +  77r  +  35r2  +  *L  ( 84  -  420r)] 

A 

m        =  108  +  252r  +  140r2  -  1008  ^5- 
id  Az 

m    =  A [-26  -  63r  -  35r2  +  ^  ( 84  -  420r)]  (20) 

14  AZ 

2 

m22  =  a2[8  +  14r  +  7l*2  +  ^2  (112  +  140r  +  280r2)] 

A 

77?  0«  =  A[26  +  63r  +  35r2  +  %■   (-84  +  420r)] 

2 

m24    =  A2 [-6  -  14r  -  7r2  +  ^  ( "28  "  140r  +  140r2)] 

2 
m,„  =  312  +  588r  +  280r2  +  1008  %■ 

JJ  A 

2    R2 
m„,  =  A [-44  -  77r  -  35r"  +  ^  ( -84  +  420r)j 

2 

m        =  A2[8  +  14r  +  7r2  +  \  ( 112  +  140r  +  280r2)] 
44  ^ 


We  note  that  if  shear  deformation  is  discarded  r  =  0,  and  if 
rotary  inertia  effects  are  discarded  R  =  0  in  Eq .  (20),  and 
m    is  greatly  simplified. 
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A3.   Initial  Stress  Stiffness  Matrix  Due  to  Axial 
Tension  Load  P 

The  following  stiffness  matrix  due  to  the  axial  tension 

load  P  is  employed  (see  2,  p.  262). 


36 

3A 

-36 

3A 

k  = 
~P 

P 
30A 

3A 

4A2 

-3A 

-  A 

-36 

-3A 

36 

-3A 

3A 

-  A2 

-3A 

4A 

(21) 
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APPENDIX  B:   PIECEWISE  CUBIC  APPROXIMATION 

Rectangular  coordinates  (x,  y)  lie  in  a  vertical  plane, 

with  the  x-axis  coinciding  with  the  axis  of  the  undeformed  barrel 

(Fig.  1).   The  axes  (x,  y)  are  attached  to  the  barrel,  so  that 

they  recoil  with  it ,  but  they  do  not  rotate  or  move  in  the  y- 

direction.   The  origin,  x  =  y  =  0,  lies  at  the  breech.   The 

barrel  is  divided  into  intervals  by  points  x,  =0,  x2  ,  x^ , 

x  ,, ,  where  x  ,n  is  the  coordinate  of  the  muzzle.   These  points 
n+1         n+1  r 

need  not  be  equally  spaced.   The  j-th  interval  is  (x.,  x.+-,). 
In  the  j-th  interval,  we  use  the  cubic  approximation, 

y  =  aj^  +  a^x  +  a^x   +  a^x  ;  x   <  x  <  x  (1) 


J  =  AJ  y   +  BJ  y     +  CJ  0  .  +  DJ  0.^,  (2) 

a     a   J     a   J+l     a   j     a   j+1 


where  0.  =  y ' ,  in  which  the  prime  denotes  the  derivative  with 
J     J 

respect  to  x .   We  set  A.  =  x.^-,  -  x..   By  algebra 
v  J     J+l     3 

A0  ■  *-3*j*l(XJ  +  l-3V    BJ0  -  ijM(3XJ  +  l-XJ) 


co--*;Sx?+i'        do  =  -^Mxj+i 


A?  =  6A.3x.  x.j.,,  B^  =  -6A,3x,  x.,-.  (3) 

1      3      3      J+l  1       J   J   J+1 


ci  "  4fxj+i(2!Ijtxjn)-  Di  "  Aj2xj(2xj+i  +  xj) 


A2  "  -3Aj3<xj  +xj+i>'   B2  "  3Aj3(xj  +xJ+i> 
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C2  =  -Aj2<*j  +2xj+l».     D2  =  -Aj2<xj+1  +  2xj> 


A3  =  2A j3 ■   B3  "  "2A j3 •    C3  =  D3  =  A j2 


AJ ,  BJ ,  CJ ,  DJ  are  constants  for  the  j-th  interval.   It  is  noted 
a   a   a   a  

that  we  get  BJ  from  A1-1  and  DJ  from  CJ  by  interchanging  x  .  ,  ,  and 
x..   This  property  simplifies  programming.   Equations  (1),  (2), 
and  (3)  define  a  piecewise  cubic  polynomial  that  is  continuous 
with  its  first  derivative. 
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APPENDIX  C: 
ACCELERATION  EFFECTS  OF  PROJECTILE 

The  following  kinematic  analysis  of  the  projectile  accele- 
ration entails  no  approximations.   Let  coordinate  axes  (x,  y) 
be  fixed  (Fig.  CI).   The  equation  of  the  tube  at  time  t  is 
x  =  x(s,t),  y  =  y(s,t),  where  s  is  arc  length  along  the  tube. 
Hence , 

2     2    n 
x  +  y   =1 
s    s  (1) 

x   =  cosG ,   y   =  sin9 
s  s 

The  location  of  the  particle  (projectile)  mass  m  at  time  t  is 

x  =  xU(t),t]  ,  y  =  yU(t),t]  (2) 

where  E,  ( t )  is  the  arc  length  s  to  the  projectile.   Hence,  by 
Eqs.  (1)  and  (2), 

sine  =  ysU(t),tJ,  cose  =  xsU(t),t]  (3) 

where  9  is  the  angle  of  the  tangent  at  point  m  CFig.  CI). 
The  velocity  components  of  the  particle  m  are 


vx  =  dt   5*s(?,t)+xt(£,t)  =  I   cose  +  xt(5,t) 
vy  =  ll  =  5ysC5,t)+ytU,t)  =  i   sine  +  ytCC,t) 


(4) 


The  terms  x    (E,  ,t)    and  yt(C,t)  are  the  velocity  components  of 
the  contiguous  point  of  the  tube.   The  normal  velocity  component 
of  m  is 
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v   =  -v   sin9  +  v   cos8  (5) 

n      x         y 

Hence,  by  Eqs .  (4)  and  (5), 

vn  =  -xtC5,t)sin9  +  yt(C,t)cos9  (6) 

This  result  shows  that  v   for  mass  m  is  the  same  as  the  normal 

n 

velocity  of  the  contiguous  point  of  the  tube.   The  tangential 
velocity  of  mass  m  is 


v   =  v   cos9  +  v   sine 
s    x         y 

=  k    +  xt(C,t)cos9  -f  yt(5,t)sin9  (7) 

or  vs  =  5  +  xtxs  +  ytys  (7a) 


where  x,x   +  y  y   denotes  the  tangential  velocity  of  the  con- 
tiguous  point  of  the  tube. 

Now  the  acceleration  components  of  the  mass  m  are 

d  x       _  d  y 


where  x  =  f[£(.t),t].   To  derive  expressions  for  a  ,  a  ,  we  must 

x   y 

more  generally  consider  the  form  x  =  f[p(t),  q(t)J.   Then, 


dx    3f  •    3f  • 
dt  =  IF  P   "3q  q 


d2x    3f  »    3f  ..    -2  32f    _•■  32f    -2  32f 

— 2        T^p+Tw7q   +   P     — ?    +   2pq  TTTTn  +  q     — 9 
dtZ         ap  3q  3p2  3p3q  3q2 


and    since   p(t)    =    £(t),    and   q(t)    =    t,    we    have 
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ax  =  ^xs(C't)+^2xss(C't)+2^xstU't)+xtt(C't)         (8) 


Similarly 


a  =  Cys(?,t)+C2yss(C,t)+2iyst(?,t)+ytt(C,t)         (9) 


The  normal  component  of  acceleration  of  m  is 


a   =  -a  sine  +  a  cos6  =  -a  y  (  5  ,  t )  +  a  x  (  ?  ,  t ) 
n     x        y         x  s   '       ysv^'y 


For  brevity,  we  may  write  x  (£,t)  =  x  ,  etc.   Hence,  by  Eqs 


(8)  and  (9)  , 


an  =  ^(xsySs  "  "saV  +  2^(xsyst  "  ysxst} 

+  (xsytt  "  ysXtt>  (1Q) 


Note  that  the  factor  E,    cancels  out  of  Eq .  (10). 
Recalling  that  the  curvature  of  the  tube  is 


1,2      2  ,1/2 
—  =  (x    +  y   )  ' 

r     ss    ss 


since  x    <<  y     1/r  «  y   .   Also,  x   s»  1.   Hence  the  term 
ss     Jss'   '     Jss        '   s 

xsyss  ~  xssys  in  Eq"  ^10)  is  aPProximately  1/r- 
The  tangential  component  of  acceleration  is 


a   =  a  cos9  +  a  sin9  =  a  x   +  a  y 
s    x        y        x  s    y  s 


Consequently  with  Eqs.  (8)  and  (9),  we  find 

as  =  ?(xs  +  ys}  +  ^(xsxss  +  ysyss} 
+2^(xsxst  +  ysyst}  +  xsxtt  +  ysytt 
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2     2 
Now,  by  geometry,  x   +  yg  =  1 .   Hence,  xgxst  +  ysyst  =  0  and 

xsxss  +  ysyss  =  °'   Therefore, 

as  =  *  +  xsxtt  +  ysytt  =  *  +  xttcose  +  yttsine       (11) 


If  we  employ  linearizing  approximations  cos6  «  1,  sin6  = 
9  =  y  .   Also,  then,  x  =  s,  x_  =  1,  xo_  =  0,  x   =  0.   Then,  by 
Eqs.  (6)  and  (7), 

vn  =  yt(?,t),  vs  =  k   +  yx(5,t)  yt(5,t)  (12) 

and  by  Eqs.  (9)  and  (10) 

ay  =  iyx(5,t)  +  52yxx(C,t)  +  2Cyxt(C,t)  +  7tt(€,t)   (13) 

a   =  I 
x 


and 


an  =  ^2yxx(?'t}  +  2^yxtC5jt)  +  ytt(?'t)  (14) 


In  the  differential  equations  of  motion  of  the  barrel,  the 

question  arises  as  to  whether  we  should  use  a  or  a   for  the 

H  y     n 

force  exerted  on  the  barrel  by  the  projectile.   This  question 
may  be  resolved  by  considering  the  tube.   When  the  tube  is 
bent,  it  is  effectively  a  curved  beam.   For  a  curved  beam,  it 
is  natural  to  deal  with  tractions  or  normal  cross  sections  of 
the  beam.   For  a  curved  beam,  we  have  the  familiar  equation 


S  +  Q  =  °  (i5) 
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where  M  is  the  bending  moment  on  a  cross  section  and  Q  is  the 
shearing  force  on  the  cross  section.   If  we  deal  with  the  large 
deflections  of  a  straight  elastic  beam,  we  have 

M  =  —  =  EIk  (16) 

r  v  ' 

where  <   -    1/r  is  the  curvature  of  the  neutral  axis.   Hence, 

o  -  -«  £ 

The  curvature  k  is  continuous  at  the  point  s  =  E, ,    but  the 
derivative  dic/ds  is  discontinuous.   The  discontinuity  is  equal 

to  ma  .   If  we  write  y  as  a  function  of  x 

n  ' 

yxx  .,3   2       , 

(i  ♦  y2)3/2  "  7=(    ~*7»      ■■ 

Hence,  with  linear  approximations,  k  =  y   .   Also,  ds  = 

^  '  Jxx 

2i/9  12 

(1    +   y    )    '  "    dx   =    ( 1    +  Tf  y    )    dx    a  dx,    and    therefore, 

Q  =    -EI    y 

Thus,  if  we  treat  the  problem  by  means  of  the  differential 
equations,  we  must  use 


AQ  =  ma   at  x  =  £ 

n 


or 


-EIAy    =  ma    at  x  =  E, 
J  xx  n 


Hence,  it  follows  from  this  argument  that  we  should  use  a  ,  as 
given  by  Eq .  (14),  to  determine  the  force  exerted  on  the  barrel 

•  9 

by  the  projectile.   Effectively  the  terms  £~  y     2£  y  t ,  y 
contribute  to  the  terms  Kx ,  Cx ,  and  Mx ,  respectively. 
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Figure  CI.       Projectile  Coordinates 
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